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Abstract 

The arguments in favor of the unified formation mechanism for both slow (Lynden- 
Bell's) bars and common fast bars are given. This mechanism consists in a certain 
instability that is akin to the well-known radial orbit instability; it is caused by the 
mutual attraction and alignment of axes of precessing star orbits (up to now, such 
a way of formation was considered only for slow bars). 

The general theory of the low-frequency modes of a disk consisting of precessing 
orbits (at different angular velocities) is presented. The problem of determining 
these modes is reduced to the integral equations of a rather simple structure. The 
characteristic pattern speeds (O p ) of the low-frequency modes are of order of the 
mean orbit precession speeds (Cl pr ). The bar- modes also belong to this type of 
modes. The slow bars have f2 p ~ f2 pr ; for the fast bars, Q p may far exceed even 
the maximum precessing speed of disk orbits (however, f2 p remains to be of order 
of these precessing speeds). The possibility of such an excess of Q p over f2™ ax is 
connected with the effect of "repelling" orbits that tend to move in the direction 
opposite to that they are being pushed. 

The preliminary analysis of the orbit precession patterns for a number of typical 
potentials is given. It is noted that the maximum radius of the "attracting" circular 
orbits (r c ) may be used as a reasonable estimate of a bar length. 



1 Introduction 



It is commonly supposed that a galactic bar can belong to one of two types — either 
common fast bars or Lynden-Bell slow bars (see, e.g., Sellwood 1993, Polyachenko 1994). 
The distinction between them is drawn along a number of lines. Firstly, bars with different 
rotation velocities have considerably different sizes: while the common bars end up at the 
corotation or 4:1 resonance, the Lynden-Bell bars end up in the vicinity of the inner 
Lindblad resonance. Secondly, it is believed that these bars are produced by entirely 
different formation mechanisms. For the slow Lynden-Bell bars, the physical mechanism 
is absolutely clear — this is the mutual attraction and "sticking" together of the slowly 
precessing orbits. Note, that in application to galactic bars, the idea of axis alignment 
was first suggested by Lynden-Bell (1979). 

Yet for the fast bars the situation remains vague and debatable. The fast bars were 
long thought to be generated due to the fast rotation in much the same way as in the 
classical incompressible Maclaurin spheroids being they are strongly oblate and rapidly 
rotating. It was shown by Toomre (1981) that in reality the galactic bar-modes have little 



in common with incompressible edge modes. Accordingly, they should have another for- 
mation mechanism. Toomre (1981) offered the so-called swing amplification mechanism, 
which currently is generally accepted for the explanation of the normal SA spirals. We 
think, however, that a simple extension of the Toomre mechanism to the SB galaxies can 
hardly be done. In particular, it is difficult to suspect the existence of the running spiral 
waves, necessary for the swing amplification, which takes place outside the bar near the 
corotation. The specific attempt made by Athanassoula and Sellwood (1986) to demon- 
strate the validity of the swing amplification mechanism on some some models of stellar 
disks appears to be somewhat artificial. 

In our opinion, it is more preferably to find such a mechanism of the common bar 
formation, which is directly linked to the instability of the central part of the galactic 
disk itself. 

Contopoulos (1975, 1977) has pointed out that the properties of the families of the 
periodic orbits in the rotating bar potentials play the fundamental role in the theory 
of the galactic bars. This is especially true for the so-called x\ family of orbits, which 
is elongated along the bar inside the corotation circle. It is plausible to assume that 
these periodic and close non-periodic orbits make up the galactic bar. One should note, 
however, that the preceding is relevant to the theory of the already existing bars, and in 
the strict sense cannot be applied the the bar formation mechanism itselfQ. 

Nevertheless, from the described picture of orbits constituting the bar, it is customary 
to make a "reasonable" conclusions about the mechanisms of bar formation, including the 
possible instability at the linear stage. The fast bar angular velocity Q p is greater (some- 
times substantially greater) than the maximum precession velocity f2™ ax . Meanwhile, it 
is intuitively suggested that the angular modulation of the precessing orbit distribution 
(i.e., a figure of the bar) should rotate with the velocity of about mean precession velocity 
of the orbits. Such a conclusion would have meant the uselessness of the Lynden-Bell 
mechanism for describing the fast bar formation. So, it is generally agreed that in reality 
the growing bar forces the orbits to change their shape and tune to the strengthening 
and narrowing bar. Moreover, it is assumed that this effect can be significant even at the 
linear stage. It is justified on the example of weak bars in the framework of the linear 
theory (Sellwood 1993): the theory shows that initially circular orbits change into slightly 
elongated ovals, with orientation relative to the bar just as the orbits of the xi-family. 

However, one can adduce two important objections as to the said in the last paragraph. 

1. The circular orbits cannot be the typical sample orbits of the undisturbed disk (ex- 
cepting for a cold disk model, which is of little interest in this case) until the bar-induced 
perturbed velocities do not exceed the velocity dispersion in the original axisymmetric 
disk. It is clear that this condition should be met first at the initial stage of the instabil- 
ity development. 

2. The intuitive conclusion mentioned above (that the velocity of the wave density 
in the system of precessing orbits, Q p , cannot exceed the maximum precession velocity 
f2™ ax ) is supported by the results of computations of the unstable modes in numerous 
models of galactic disks. These computations show that in all cases the pattern speed Q p 
is smaller than the maximum star angular velocity fi max in the diskQ. But this analogy 

1 This is especially true for the linear stage of the bar- forming instability (if the bar arose due to the 
growth of the instability): the orbits considered by Contopoulos are obviously trapped by the potential 
of the bar. The trapping process is, of course, the non-linear phenomenon. 

2 Note that this property of the oscillation frequency spectrum of the gravitating disk is not yet found 



appears to be somewhat loosely after the discovering of the possible "donkey" behavior 
of the star orbits by Lynden-Bell and Kalnajs (1972), when the orbit accelerates when 
held back and slows down when urged forward. Such a behavior takes place under the 
Lynden-Bell condition, dQ pr /dL\j < (Lynden-Bell 1979), where L is the star angular 

momentum, Jf = J r + L/2 is the Lynden-Bell adiabatic invariant , J r is the radial action. 
As it is shown below, if the orbits with a "donkey" behavior play an active role in the 
bar-forming instability, the bar angular velocity can exceed f2™ ax . 

The last argument allows to treat the fast bar-modes as the corresponding density 
waves of the precessing orbits, in the full correspondence with the theory of the Lynden- 
Bell bars. Note that Kalnajs (1973) was first who called attention to the possibility of the 
freely precessing orbit alignment (and thus the bar formation). The important point in 
his theory of kinematic waves in the case of near-circular orbits was the independence of 
the precession velocity £l pr (r) = Q(r) — k{t)/2 (fi( r ) from radius (here n(r) is the epicyclic 
frequency, k 2 = AQ 2 + rdQ 2 /dr). Note that the same condition Q pr (r) m const is required 
in the Lynden-Bell (1979) more general consideration of the problem. 

In reality different stars precessing with different velocities. But only important point 
is that the angular precession velocities (and also the pattern speeds of bars (Q p ) in the 
majority of cases) are substantially smaller than the typical star azimuthal frequencies Vl 2 
and especially the radial frequencies Qj.. If for some orbit the inequality 

\n p - JVIM < i (i) 

holds, this whole orbit (not individual stars) can be considered as an elementary object in 
the interaction with the bar gravitational field. If the condition (Jl|) holds for the majority 
of disk orbits, participating in the bar formation, we can consider the processes (e.g., bar- 
instability) in the model disk, consisting of the set of the precessing stars. As it is noted 
by Lynden-Bell (1979), the condition (|]) means the conservation of the adiabatic invariant 
Jf = J r + L/2, that reduces the problem of the bar- instability to the one dimensional 
problem: one have to follow only the variation of the azimuthal positions of the orbit 
major axes under the action of the gravitational attraction from the bar. 

Actually, we solve the problem of bar density wave of orbits with different precession 
velocities, that is quite analogous to more familiar problem of the star density waves 
(e.g., of the bar-like shape) in the differentially rotating disks. Seemingly, not all the orbits 
participating in the bar-mode formation, meet the condition ([]]) with margin, especially for 
the fast bars. But even for the latter usually |fi p — Cl pr \ ~ Q pr ~ uq, where u G ~ ^GMd/a? 
is the characteristic gravitational (Jeans) frequency (G is the gravitational constant, 
is the mass of the active disk, a is its radius). Then when Q pr /Qi < 1 we are still within 
the bounds of the condition (p. Even in case of more weak inequality Vt pr /Vli < 1 (e.g., 
several times smaller but not at some orders) the suggested model will likely provide the 
correct qualitative answer. This is at any case not worse than, for example, the analysis 
of the spiral structure of the galaxy M 33 by Shu et al. (1971) by using the WKB formulae 
from the well-known theory of Lin and Shu, applicable, strictly speaking, to the tightly 
wound multi-turn spirals. 

A few words about the content of the paper. In the Section 2, we descibe in the general 
form a model of the precessing orbits and derive the basic equations for the model, useful 



for the general case. 



for the analysis of the disk low frequency modes of interest for us. The use of the action- 
angle variables Ji, J 2 = L and Wi, w 2 is the easiest way for derivation of these equations 
from the general kinetic equation. The substantial simplification appears with introducing 
a slow angular variable w 2 = w^ — Wxjl and averaging over the fast angle variable w\ (i.e., 
over the radial stellar oscillations). In the simplest cases, the problem can be reduced to 
the analysis of a rather simple dispersion relations. In some more general cases, we obtain 
the integral equations of different degrees of complexity. But even in the most general 
form, the obtained integral equations for this model is much simpler that the immense 
integral equations for the normal modes of the stellar disk, obtained by Kalnajs (1965) 
and Shu (1970). Remind that the use of these general integral equations (not counting the 
iV-body methods) was the only possibility for analysis of large-scale modes, as opposed 
to incomparably simpler problem of tightly-wound spirals . This is especially true in 
regard to the bar-mode. As discussed above, the real simplification has come through 
the analysis of low-frequency modes. Note that the resulting description of gravitating 
systems is similar to the drift approximation in the plasma physics (see, e.g., Chew et. al., 
1956), but for orbits of essentially more general type. In our opinion, the most important 
advantage of our approach consists in the fact that it makes clear the simple physical 
mechanisms of the instability processes developing in a disk. Unfortunately, a possibility 
of revealing these physical mechanisms under the use of the general integral equations 
by Kalnajs or Shu would be practically impossible, and the same is true for the iV-body 
simulations. 

In Sec. 3, we analyze the dispersion relation for a model disk consisting of orbits of 
two different types that differ by their precessing speeds (O^ and fi^; Sv® > and, 
generally speaking, also by a sign of the Lynden-Bell derivative (dQ pr /dL) Jf = Q! pr .For 
the case when the derivatives (Q! )^ and (Vt' pr )^ have opposite signs (i.e., a disk contains 
both "attracting" and "repelling" orbits), the pattern speed of the unstable modes KeQ p 
may be more than (f^,.)^, i.e. the maximum orbit precessing speed in a disk. We show 
that at the same time VL p l ) < Ref2 p < flW when (fl' )W and (f2 pr )( 2 ) have the same sign 
(positive for instability). We consider this result as the important argument in favor of 
the unified formation mechanism both for slow and fast bars. 

In Sec. 4, the results of computation of the precessing speeds Q pr and the most 
important for the theory quantity, the Lynden-Bell derivative Q' pr , are given for a number 
of typical potentials. In particular, the regions on the Lynden-Bell Jf,L plane where 
Q' pr > and Q' pr < are found. The natural suggestion is made that a bar forms by 
"attracting" orbits with Q' pr > 00. Then the bar length (more exactly, the radius of 
the bar's end) should be equal to the maximum of the apogee radii (r max ), among the 
sufficiently occupied orbits from the region where fl' pr > 0. So, for calculation of lb, one 
should know not only the general pattern of orbit precessions determined by the potential 
$o(?") but a specific equilibrium distribution function as well. As a first approximation, 
one can take the estimate 1^ ~ r c , where r c is the radius of the circular orbit at which 
Q' pr = (Q' pr > for r < r c , and Q' pr < for r > r c , at the circular orbits). These radii 
r c are computed for all the potentials considered. Note that our determination of the 
bar length has nothing to do with common determinations that link this length with a 
location of one of the resonances (CR, ILR or 4:1). It is clear that our bar length may 

3 Note that in the formation process for all the mode, i.e. both the bar and adjacent spirals, and not 
just the central bar, the "repelling" orbits with il' pr < may also take part. Moreover, this may be 
important to explain the fast bar phenomenon (as discussed above). 



take a great variety of values depending on the specific potential and distribution function 
(accidentally, they may fall near one of resonances). 

In Sec. 5, we shortly formulate the most important conclusions and discuss some 
immediate prospects for a work in this field. 



2 Basic equations for the low-frequency modes 

For derivation of the basic equations of the theory, it is most conveniently to use the 
action-angle variables I = (ii, J 2 ) and w = (wi,w 2 ), which suitably takes account of the 
double periodicity of stellar motion in the equilibrium potential. Note that Ii = I r (J r is 
the radial action), I 2 = L (L is the angular momentum). We start out with the linearized 
kinetic equation in its usual form (see, e.g., Fridman and Polyachenko 1984): 

+ (2) 
dt dw% dw 2 dli dwi dl 2 dw 2 

where /o(I) and /i(I, w, t) are the unperturbed and perturbed distribution functions, $1 is 
the perturbation of the gravitational potential, and Q 2 are the frequencies of the radial 
and azimuthal oscillations of stars in the equilibrium potential $o( r )> — dE(T)/dIi (E is 
the energy in terms of I, i = 1, 2). The change of variables w 2 = w 2 — Wi/2, W\ = W\ in 
(1) yields 

df M . n ,, n df (df ldf \ 

— — h im\l vr t + ih— — = — — — h Jmf — — — , (3) 

dt pJ dw x dhdwi \dl 2 2dh)' y! 

where we have assumed that the perturbations are proportional to exp(imw 2 ): 

$1 = Qexp(imw 2 ), fi = f exp(imw 2 ), 

m is the azimuthal index (an even integer), and Q pr (E, L) = Q 2 — is the precessing 
speed of an orbit with energy E and angular momentum L. If we also transform from the 
action (/]_, J 2 ) to (E, L) (the equilibrium distribution function is usually given just in the 
variables E, L), we obtain the linearized kinetic equation in the form 



OF . _ 8F dF 9$ . ^ fdF dF 



where F (E,L) = F (ii,I 2 ). Note that in this form, Eq. (^) also holds for a quasi- 
Coulomb potential i- e - ; one that is due principally to a large central mass. All that is 
necessary then is to redefine w 2 and Q pr : w 2 = w 2 — Wi, Q pr = Q 2 — f^. 

Then, neglecting self-gravitation of the system and the quadrupole moment of the 
central object, we have fl pr = 0, as it should be for closed Keplerian orbits. Here the 
azimuthal index m can be either odd or even. 

As it is shown by Lynden-Bell (1979) (see details in the Introduction) the most ad- 
equate variables that should be used instead (/j, 7 2 ) or (E, L) are Jf — I\ + J 2 /2 and 
L = I 2 . Changing in (^) to these variables, we obtain the kinetic equation in the most 
convenient form for studying the low-frequency modes: 



(5) 



where J-o{Jf, L) = fo(h, h) and taken into account that dj^jdli — ^df /dli = dF /dL. 

We assume that the rms deviation of the precession rates about the mean Q pr , given 
by Att pr = [(Qpr — VLpr) 2 ] 1 / 2 and the typical gravitational frequency ujq are both small, 
Cl pr , loq <C Qi. Note that ujq is of the order of the mean Jeans frequency of the system: 



ojq ~ yGM d /a 3 , where M d and a is the mass and the radius of the active disk, respec- 
tively The above inequalities can be justified, e.g., if one assume that we are dealing with 
a system of stars, that have close precession speeds, within a massive halo which, while 
remaining unpertubed itself, furnishes the dominant contribution to the equilibrium po- 
tential $ - However, it should be remembered that the role of halo should not be played 
only by the real massive spherical component of the galaxy. Indeed, since in general the 
extent to which different groups of orbits are involved into the instability process can 
differ greatly, one can consider as the first approximation that the active group of stars is 
immersed in the massive halo of other disk stars. Under these circumstances, then, there 
may exist a low-frequency mode (oc exp(— iujt), with uj = uj — mQ pr ~ ujq, AQ pr ) in a 
coordinate system rotating at angular velocity fl pr such that the slow precessional disper- 
sal of orbits is canceled by their mutual gravitational attraction. It would be natural to 
suppose that if self-gravitation were to win out over the dispersion in orbital precession 
speeds, an instability should develop that could eventually deform the system (under the 
influence of the largest-scale growing modes). It is clear, however, that even in a system 
with essentially radial orbits, this holds true only if the torque that alters the orbital 
angular momentum of the stars forces their orbital precession speeds to change in the 
same direction (see the Introduction and, e.g., the formula (|19]) below). 

We use the perturbation theory to derive the desired solution for the low-frequency 
modes. Let T = .F^ 1 ) +J r( - 2 \ where J 7 ^' corresponds to the permutational mode obtained 
from ([5]) by neglecting terms proportional to AQ pr and $ oc G: uj = (or uj = mCl pr 
when VL pr ^ 0), dJ-'^/dw 1 = 0, i.e. J 7 ^ = JF^^Jj, L) is an as-yet arbitrary function of 
the integrals of motion, which we will subsequently specify by requiring that the solution 
of the next approximation be periodic. The equation for T^ 1 takes the form 

- imfW + zmlV^ 1 ) + = + ( 6 ) 

owi oJf owi oL 

Averaging (^) over Wi, from to 2it, and bearing in mind the periodicity of the functions 
JF( 2 ) and $, we have 

f)T 1 2 r 

-{uj- m5n pr )^ w / $d Wl , (7) 



dL 2tt 



o 



where 5il pr = il pr — Q pr . 



Invoking the Poisson equation, some minor manipulations yield 

$ = -G J dJ'^(J') j dwT(r,r' , ip' - varphi) exp[im(w' 2 - w 2 )\, (8) 
where dJ' = dJ'jdL', dw' = dw[dw'2, T is the Green function: 



T = — , ri2 = [r 2 + r' 2 — 2rr' cos(<// — f)] 1 ^ 2 - 

r 12 



(9) 



Equation (H) is an integral equation for the function <E>(J, wi),) if one takes into account 
the expression (|7|) for J 7 ^ through $. The coordinates of stars r, ip, r', ip 1 in (||) and (||) 
must be expressed in terms of J, J', w, w', where 

r = r(J, w), r' = r'(J, u^), w 2 — w 2 = (w' 2 — W2) — (w[ — wi)/2, 

ip' — ip = Sip = w' 2 — w 2 + 0(J, J', Wi, w[) 

and we refrain from writing out the expression for 0. Since the perturbed potential can 
always be written out as 

$i(r, (p) = 5 ) i(r) exp(irrnp) = <£>exp(OT«D 2 ), 

we have 

$(r, ip) = $i(r) exp[zm5(J, Wi)], 

(5 = ip — W2 is a known function of J and Wi). In actual fact, then, ([8]) is an integral 
equation for the unknown function $i(r) of only one variable. This equation can be 
rewritten in more symmetric form. The right-hand side of (H) depends on W\ only through 
T and exp[im(w 2 — 102)]- So, by averaging (H) over wi, we obtain for the function 

1 2?T 

X(J) = $ = — / 



the following integral equation 

x(J) = ^ / dJH(J, JO dM 2 /9 n<\ xW (10) 

where 

il(J,J') = y dw 1 (iw / 1 (i5w2r(r, r', 6(p) exp(zm5u7 2 ) exp[— im{w' l — w{)/2\, (11) 
and <5w 2 = w' 2 — w 2 - 

Physically, n(J, J') is proportional to the torque 5M acting upon some selected (test) 
orbit with the action J resulting from all orbits with fixed action J'; these all have the 
same shape, but their major axes are oriented in all possible directions: 

8M oc -imGexp(imw 2 )Tl(J, J')J r(1) (J')^J'- ( 12 ) 

For a quasi-Coulomb field $o( r ) ; instead of the Lynden-Bell integral Jf = Ji + 1 2 /2 
we have to use Jf = I\ + 12 in ( |10|) and instead of exp[— im(w[ —wi)/2] (with m required 
to be even) in (|TT| ) we have to use exp[— im(w' l — wi)] (with arbitrary m). 

One might hope to reduce ( |10D to one-dimensional integral equations in two limiting 
cases: 1) when the distribution function .Fq(J) is close to a delta function in L near some 
value L ; we will be commenting on this circumstance, writing JF = A±(Jf,L — L ) rs 
8(L — L )ip (Jf); 2) for the systems with near-circular orbits, whereupon f = A 2 (/i, I2), 
where A 2 ~ 5(Ii){po(l2). Below we will restrict ourselves mainly with the first case. The 
second case is technically somewhat more complicated; we will study this important case 
elsewhere. 



Thus we assume Tq = Ai(J/,L — Lq). Using the fact that II and x vary in L' only 
a little over the characteristic scale length of the function (dJ-'o/ 'dL) /[u — m5Q pr (J'f, L')], 
we can reduce Eq . C[TQ|) to an integral equation for the function of one variable if>(J) = 
x(J, L' = Lq); for brevity, hereafter we omit the index "/" in the Lynden-Bell integral Jf. 

1>{J) = ^J dJ'P(J, J')So(J'MJ'), (13) 

where 

P(J, J') = n( J, J', L = Lo, U = L ), (14) 

5u(J) _ f dL > 9W,L')/dL> 
J to — m50 pr ( J', U) 

More convenient form for the function Sq(J') is obtained after integration (|T4| ) by 
parts 

qn/ , f Jr/ W,L')dQ pr (J',L')/dL> 

S (J) = -mJdL [{D _ m6nAJfjLI)]2 , (16) 

where the derivative dQ pr /dL' can be taken at L' = L . 

If we are dealing with near radial orbits, we can put L = when we calculate 
P(J, J'). Furthermore, 5ip = ip' — tp ~ wf, — u>2 for such orbits, so the function II can then 
be substantially simplified: 

II(J, J') = J tftM«^m[r(J,wi)MJ>'i)], ( 17 ) 

where 

2tt 

J m (r,r') = — I daT(r, r', a) cos ma. (1£ 



: — / daT(r,r', a) cos ma 
2ix J 



2tt 
o 

When the orbits are exactly radial ("cold" system), i.e., Tq = S(L)<p (J), 

Tfi 

So(J') = --A(J')MJ'), (19) 



where 

dQ pr (J',L') 



A(J') 



dU 



(20) 

L'=0 



Accordingly, the integral equation (0), which then describes the radial orbit instability 
in a cold system, is 

CtrrP' f 

*KJ) = "2^5 J dJ>P (J, J')A(J')MJ'W)- (21) 

It can easily be shown that J m (r,r') defined by QTSp is a positive function. Conse- 
quently, so is II(J,J / ) from (PI) , and most importantly, so is P(J,J') from (fUj) when 
L = 0; just the function P(J,J') enters into the derived equation (|2~I|). This can be 
made explicitly by expanding the function (r 2 + r' 2 — 2rr' cos a) -1 / 2 in a harmonic series, 
we obtain 

J m {r,r') = VF„(r,r')- / daP n (cosa) cos', 



ima, 



where 

F n (r,r') = r^/r™ +1 , r < = min(r, r'), r > = max(r, r'), (22) 

and the P n are Legendre polynomials. Now the positive defmiteness of J m follows from 
the fact that P n (cosa) can in turn be expanded in cosines of multiples of the angle, with 
positive coefficients (Gradshtein, Ryzhik 1965): 



x (2n-l)!! 

P„(cosa;) = : — ; — 

v ; 2"-%! 



with all > (a prime indicates that the parity of k and n must be the same), and 
from 

1 f j Omn 

— I da cos ma cos na 



In 

cosna H cosm — 2)a + 

1 2n - 1 V ; 



= Afr , cos ka 

k 



TV J 2 


We finally obtain a convenient representation for J m (r, r') in simple series form: 



4(r,r')^E'F n (r,r')4")>0. 

n>m 

Given the positivity of P(J, J'), the sign of the integrand in (0), and therefore the 
sign of u 2 (i.e., the stability or instability of a system with purely radial orbits), will 
depend on the sign of A (J') as defined by (|T9|) . If A > for all orbits in the system under 
consideration (in other words, for all values of J or what amounts to he same thing in the 
present case, for any value of the energy E of radial stellar oscillations), then u 2 < 0. As 
a result, radial orbits are unstable with A > 0. On the other hand, provided that A < 0, 
it is the purely oscillatory mode. The most compact formula for computing A(E) is 

dx 1 




*o\J x 2 [l - %{x) / E} 1 / 2 r 

A{E) = j2Eyf2 r \ j x 







[2E - 2$ (o;)] 1 / 



2 



The inequality (dQ pr /dL)\L =Q > is merely a necessary (and in no way sufficient) 
condition for the radial orbit instability, and in particular, for the formation of a bar. 
The insufficiency of this criterion is immediately obvious from the fact that the retarding 
torque due to the bar can turn out to be ineffectual in the face of large orbital precession 
speeds. To derive valid conditions for bar formation, it is necessary to solve the problem of 
stabilization of the radial orbit instability by some finite dispersion of the orbit precession 
speeds [having demonstrated once again the predominance of the bar mode (m = 2)]. 
The bar formation criterion is in fact none other than the condition for the bar-mode 
instability of the type under consideration. We therefore now proceed to derive the 
stabilization conditions for systems with near-radial orbits. 

For defmiteness, we assume the Maxwellian distribution in L 

/o = -Tjrr e M-L 2 /L 2 T )MJ), (23) 



where Lt is the thermal spread. If we then assume in (^) that u> = uj = 0, L = 0, 
5Q pr = Q pr ~ A(J)L, we will then have for the system's stability boundary 

So(J') = 2<f (J')/mLlA{J') 

so that the integral equation ([TJ) becomes 

m = ^jj dfpu (24) 

This equation is almost the same as Eq. (PI). A comparison of these two equations 
suggests a simple relationship between the instability growth rate 7 for a system with 
purely radial orbits, 7 2 = — u 2 , and the minimum dispersion in orbital angular momentum 
required to supress that instability: 

(L T ) min = 2 l ' 2 1 /mA, (25) 

where A is some mean over the stellar orbits with different energies E. 

The relation ( 25|) acquires a precisely defined meaning when all stars have almost 



the same energy, E Eq, since we can then take A = A(Eq). If in fl2"B| ) we go to the 
distribution over precession speeds, Q pr = AL, we then obtain a more obvious relation in 
place of (p5f): 

(n pr ) T = 2 1 / 2 7 /m, (26) 

where (fl pr )T denotes the thermal spread in precession speeds, and the growth rate 7 is 
given in the form 7(771) to emphasize that in general it depends on the azimuthal index 
m. Since 7(771) is only a weak function of 7770, it follows from (^) that the most difficult 
modes to stabilize (and in that sense, the most unstable) are those with the smallest 
possible 777. For almost radial orbits, we have m m \ n = 2, which corresponds precisely to 
formation of an elliptical bar out of an initially circular disk. All the modes with odd m, 
particularly the m — 1 mode, are suppressed in this case, as two oppositely directed (but 
equal) moments of forces would act on the two halves of an elongated orbit. The forces 
break, but do not rotate such "needle" orbits. 

For near-circular orbits in a potential close to that produced by a central point mass, 
however, the m = 1 mode is immediately become dominating. 

In conclusion to this Section, we give the integral equation of a type (fl3|) in a form, 
suitable for computation of the eigen frequencies of the low-frequency modes of a stellar 
disk with the equilibrium distribution function f (E,L), provided that a disk contains 
substantial fraction of elongated orbits: 

■S'max 

f(E 1 )= J K(E 1 ,E 2 )f(E 2 )dE 2 , (27) 

-^min 

where the kernel is equal to 

imax . , , a a 

K(E 1 ,E 2 ) = --^— I *» I *' * > OrfaLjdZ* J Jdxdyp^\x)p^\y)J m (x,y), 

l{ 1) Q [Hp ilpr ) Q Q 

(28) 

4 For example, we have 7(777) oc m 1 / 2 for m ^> 1. Then, the equation (|l7| ) is simplified since J m « 2<5(r— 
r')/m for m 3> 1. Thus, one of the integrations in ( fl7| ) can be performed. The asymptotic expression for 
J m is most conveniently derived directly from Poisson's equation, assuming that rn 2 $i/r 2 ^> \d 2 $i/dr 2 \, 
\r~ 1 d$i/dr\ in the latter. 



Such a form of the integral equation suggests that the torque of forces produced by 
attraction of two elongated orbits can be approximated by considering each real oval 
precessing orbit as a "needle" , which coincides with a major oval axis; the linear density 
of a needle is p\ E ^ = l/v r (E) = 1/ \J2E — 2& (r), E is the star energy, v r is the star radial 

velocity, M 1 (E) = J p\ [r)dr is the half of needle mass, 2a is the needle length, and 
o 



n'(E,L) 



oa 



pr 



■If 



~dT + pr ~dE~- 



By using the equation ( p7|) Polyachenko (1992) had computed the "abnormally" low- 
frequency bar-modes, which were found by Athanassoula and Sellwood (1986) in their 
iV-body study of the linear stability of some exact phase models of stellar disks. Those 
frequencies were "abnormally" low compared to the frequencies of "standard" fast bars, 
obtained by them for the majority of the models studied. Actually, the pattern speeds 
of the low-frequency modes are approximately equal to mean orbit precession speeds in 
a central disk region. Hence it follows that the instability of elongated orbits occurs. 
Figs. |I|a, a' show a typical orbit participating in the slow bar- mode instability; this orbit 
corresponds to mean values of energy and angular momentum of stars over the region of 
the bar location. As one can see, the orbit is strongly elongated; this fact justifies the use 
of the equation (|27|) with the kernel (P8|). The computed growth rates (Polyachenko 1992) 
are in good agreement with those obtained in the paper by Athanassoula and Sellwood 
(1986). On the other side, Figs, figlb, b' show the analogous typical orbit for the models 
in which only the fast bar mode developed. This orbit is substantially more round, than 
the orbit in Figs. |l]a, a'. So, for obtaining the low-frequency eigen modes in such models, 
the use of the general integral equations (||) or ( |T0"D seems to be more adequate. However, 
we think that even the integral equation (|27|) , being certainly a rough approximation 
for these models, can provide a satisfactory numerical agreement with iV-body results of 
Athanassoula and Sellwood (1986). We plan to study all these problems elsewhere. 



3 Bar-mode in the model two-component disk 



Let us represent the dispersion relation that can be derived from (|27|) for "one-component" 
system with the equilibrium distribution function f = A5(E — Eq^)5(L — Lq )> i n the 
following form: 

1 + <Ay^°- (29) 

where Qi = Q pr (E^\ Lq ) and g\ oc Q'^Eq 1 ^, L^) denotes the corresponding coefficient; 
its explicit expression can be easily obtained from ( |27|) after substituting a given 5- 
distribution function. In (p9|), the instability corresponds to g± < 0. 

Below we consider the case of the two-component system, with the distribution func- 
tion f = A5(E - E^)5(L - 4 1} ) + B5(E - E [ q ] )5(L - Lf ] ). Then some simple manip- 
ulations with the equation that is obtained after substituting this distribution function 
into (pTD lead to the dispersion relation: 

1 , 9i 92 . 9i92 „ , Qn x 





Figure 1: Typical orbits of stars in the Shuster potentials for the models studied by 
Athanassoula and Sellwood (1986): a- models with the smallest mean precession speeds 
and angular momenta; a' - the same as a but in the reference frame rotating with the 
precessing orbit; b - majority of models; b' - the same as b but in the reference frame 
rotating with the precessing orbit. 



where g 2 oc Q pr (E^\ L^) is analogous to the coefficient g 2 for the first component intro- 
duced earlier, Q 2 = ^ P r(E^\ iff), and the designations are used: 

7-2/ p (2)n 

1 A) (A) ^0 j /oi\ 



IoiE®,^) =j J dxdyp^\x)p^\y)J m (x,y), (32) 







It turns out that the possible locations of Re fl p for the unstable roots (7 = Im Q p > 0) 
essentially depends on the signs of the coefficients g\ and g 2 . 

Let us prove first of all that for positive g\ and g 2 the inequalities Oi < Refi p < Q 2 
occur. To do this, it is sufficient to calculate the imaginary part of the left side of the 
dispersion relation (23). We have 

Al ' 2 = Im (^-^ 2)2 = - (Af i2 + 7 ^ ' A ^ = Re ^-^; 



5 = Im (n p - n,nn p - n 2 y ~ " (A f + 7 T(Ai + 7 t [( 2 + 7 } 1 + ( 1 + 7 )Aa] ' 



Assuming Ref2 p > f2 2 we have Ai > 0, A 2 > 0. So in this case the imaginary part of the 
left side of the dispersion relation (23) would be negative: 

A l + A 2 + B < 0, 
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Figure 2: The trajectories of the unstable root in the two-component model with opposite 
signs of g\ < and g 2 > 0, for a number of values of the coefficient a = 0, 0.2, ... 0.7. 
The absolute value of g\ increases along the trajectories from the starting point g\ = 0. 

if one takes into account that 7 > for the unstable roots of interest and, besides, a > 
as it follows from the Cauchi-Bunyakovsky inequality (positive definiteness of the weight 
function J 2 (^, 2/) was proved in the preceding section). 

Similarly, for Ref2 p < Qi all inequalities are reversed (of course, except for a > 0): 
Ai < 0, A 2 < 0, so 

A x + A 2 + B > 0. 

A completely different type of situation occurs for the case when the signs of g\ and g 2 
are opposite, i.e. the disk contains the orbits with a "donkey" behavior; for definiteness, 
we assume that g 2 > 0, g\ < 0. Fig.[I| show the trajectories of motion of the unstable 
root on the complex plane fl p , for a fixed g 2 = 0.05 2 > (Q± = and fl 2 = 0.25 are also 
fixed) and negative g\ that vary along each trajectory from gi — to some (gi) m in, at 
which this root become stable too. Different trajectories correspond to different values of 
the parameter a. The most important fact following from these calculations is that the 
excess of Ref2 p over can be quite significant (about 1.5 times in a given example). 

4 The patterns of orbit precessions in some typical 
potentials 

Fig. [^a-|6|a show, on the Lynden-Bell plane (Jf,L), the constant value curves for the 
derivative Q' pr = (dfl pr /dL) Jf , for a number of the commonly occuring potentials $o( r )- 
In the parallel Figs.|3|b-[(]b, we give, for the same potentials, the angular velocities of stars 
at the circular orbits fi(r), tthe rotation curves Voir) = rQ(r) and the precession speeds 
of the nealy-circular orbits Q pr (r) = Q(r) — k{t)/2 (n(r) is the epicyclic frequency). 





Figure 3: a - the pattern of orbit precessions for the isochrone model at the Lynden-Bell 
(Jf, L)-plane. The strait line Jf = L/2 corresponds tocircular orbits. The curves are the 
isolines for the Lynden-Bell derivative Q' . b - the angular disk velocity fl(r), the rotation 
curve Vo(r) and the circular orbit precession speed Q pr (r) for the isochrone model. 

The first pair of these figures (Figs.^a, b) correspond to the isochrone potential 

1 



1 + VT+r^ 

This case was earlier considered by Lynden-Bell (1979). For this potential, the frequen- 
cies Qi(Jf,L), D, 2 (Jf,L) as well as the quantities of interest fl pr (Jf,L) = Q 2 (Jf,L) — 
Qi(Jf, L)/2 and ft' (Jf,L) can be obtained analytically. We used the example of the 
isochrone potential as a test model for our general scheme of computation of Q pr and Q' pr 
for the arbitrary potential $o( r )- 

Apart from the isochrone model, we carried out the computations for the Shuster 
potential 

(Figs. (|a, b), the logarithmic potential $o( r ) — h 17 ", corresponding to the flat rotation 
curve Vo = const (Figs.|5|a, b), and for the potential of the exponential disk, 

%(r) = rh(r/2)K (r/2), 




Figure 4: The same as in the Fig. ^ for the Shuster potential. 

where I\ and K are the corresponding Bessel functions (Figs.|6]a, b). Qualitatively, 
Figs.|3]a,f|a,|]a and |3|b,|4]b,|6|b are similar, but they differ greatly from Figs.^a, b for the 
logarithmic potential (and it is not unreasonable). In our opinion, the most interesting 
information that one can extract from Figs.|3|a,|]a,|6]a are the critical values of the angular 
momentum (L c ) that separate the regions of "attracting" (when L < L c ) and "repelling" 
(when L > L c ) near-circular orbits. The critical radii (r c ) corresponding to these values 
of L c are shown in ^3,^3,^3. It is natural to take r c as an estimate for a bar length 
lb that forms as a result of the bar-instability under consideration (a least at the linear 
stage). As one can see, r c is always more than the radii corresponding to maxima of the 
rotation curve Vo(r) or the function Q pr (r). Under the condition commonly used that 
a fast bar ends near the corotation, we obtain from lb ~ r c such an estimate of the bar 
pattern speed, for the Shuster potential: Q p ~ 0.23 (see Fig.|]b). The corresponding eigen 
frequency u = 2Q P w 0.46 is typical for the majority of models studied by Athanassoula 
and Sellwood (1986), just for the Shuster potential. Note that these bars are traditionally 
considered as "fast" ones, keeping in mind their non-Lynden-Bell formation mechanism. 
We see, however, that the fast bars can likely be formed by the same mechanism as for 
the slow bars. Let us point out some other interesting regularities that are common 
for the patterns of the orbit precessions in potentials $o( r ) of a type corresponding to 
Figs.|a,|a,|a: 
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Figure 5: The same as in the Fig. [$] for the logarithmic potential. 



1. The precession speed (fi^ ax ) corresponding to the maximum of the curve Q 



pr \ 



(at the circular orbits) is the absolute maximum for the precession speeds Vt pr (Jf,L) of 
arbitrary orbits; (fi™ ax ) f° r three models are given in the first line of the table. 
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2. (fi'™ ax ) are approached at the circular orbits in the disk center (Jf — > J r — > 0, 
L — > 0); the values of (f2'™ ax ) for three models are given in the third line of the table. 

3. (fi'™ m ) are approached in some points at the circular orbits; the values of (Q' 
are given in the fourth line of the table. 

The pattern of the orbit precession in the logarithmic potential (typical for the major 
parts of many spiral galaxies) is entirely different from the other cases considered (see 
Fig. ||a). The most interesting fact here is that Q' pr < only within a rather narrow sector 
of the (Jf, L)-plane, adjacent to the line of circular orbits Jf = L/2. 

The massives of the values of functions fl pr (Jf,L), Q' pr (Jf,L) obtained above will be 




Figure 6: The same as in the Fig.|3] for the potential of the exponential disk. 

used later under studying the integral equations (||), ( fL0|) for deriving the eigen bar- modes 
(first of all, the stellar models in the Shuster potential). 



5 Conclusion 

Let us formulate and discuss some conclusions from the theory above. 

1. We advanced a number of arguments in favor of universality of the mechanism 
that may be responsible for formation of both the slow and fast bars. The essence of this 
mechanism is naturally formulated on the basis of representing a stellar disk as a set of 
precessing orbits. Such a concept is adequate to the problem under consideration since 
the bar pattern speeds (including "fast bars" ) are significantly less than the characteristic 
frequencies of oscillations of individual stars (Qi and f2 2 ), but they are just of order of the 
orbit precession speeds (fl pr )- In such a disk, we seek the unstable normal modes (first 
of all the bar-mode) as the density wave of precessing orbits that runs with some speed 
Q p without any deformations despite the fact that different orbits precess with different 
speeds (analogously to differentiability of a disk rotation when the usual concept of a 
galactic disk set of individual stars is used). 

The unstable bar-mode forms if a central region of a disk (a location of a future bar) 
contains a sufficiently massive group of "attracting" orbits that satisfy the Lynden-Bell 



condition Q' pr > 0, and, besides, the precession speed dispersion of these active orbits are 
not too large (otherwise, the orbits will run away from the region of perturbation under 
the influence of the "thermal" motion). The last condition is natural for the Jeans nature 
of the instability under consideration. The exact criteria of instability can be obtained by 
solution of the basic equations derived in Sec. 2; for some simple cases, the corresponding 
dispersion relations are given in an explicit form. The pattern speed of a bar fl p depends 
significantly on the extent to which the "repelling" orbits with VL' pr < (the orbits with 
a "donkey" behavior) take part in the bar-formation process. If such orbits are hardly 
dragged in the bar-formation process, then Q p ~ Cl pr ; just such bars would naturally be 
named as the "slow" bars. But if the role of the "repelling" orbits is essential, we can 
obtain the "fast" bars with Q p that significantly exceeds O pr (just the same, Q p should be 
of order of Vt pr as before — if we want to remain within the framework of our theory). In 
Sec. 3, such a possibility is demonstrated on the simplest example of a two-component 
disk model. Thus, from the point of view under consideration, distinctions between the 
slow and fast bars are mainly quantitative, but they do not differ fundamentally from each 
other: both bars form under the action of the same physical mechanism. It is worth noting 
that the Jeans mechanism (including one under consideration) is always best natural and 
suited to the gravitational problems. 

2. So far the theory above was confirmed only by the calculations of the lowest- 
frequency modes of the disk models in the Shuster potential when the results are compared 
with the iV-body results by Athanassoula & Sellwood (1986). These modes obviously 
correspond to the slow bars: for them, Vt p ~ Q pr . The most interesting prediction of the 
theory (its validity for the fast bars) would be tested if we shall be able to prove that 
the eigen frequencies of the modes computed from our integral equations (see Sec. 2) 
coincide with the frequencies of the fast bars derived from the iV-body simulations (in 
particular, with the majority of bar-modes from the paper by Athanassoula & Sellwood 
cited above. This problem will be the subject of study in the immediate future. In the 
present paper, we restricted ourselves only to some positive facts that result from the 
general analysis of orbit precessions for a number of potentials (Sec. 4). These facts 
correlate with a possibility of the fast bar formation by alignment of "attracting" orbits. 
First of all we noted that the maximum radius of nearly-circular orbits (r c ) with Q pr > 
(r c may be taken as a natural estimate for a bar end lb) is significantly more than a size 
(~ r m ) of the solid rotation region, for all the reasonable models. Moreover, if one takes 
(as it usually does) that lb ~ tcr (where tcr is the corotation radius), then for lb ~ r c 
one can obtain the estimate uo w 0.46 for a typical eigen frequency u of bar- modes in the 
Shuster potential; this estimate is in good agreement with the corresponding results of 
Athanassoula & Sellwood (1986). 
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